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ABSTRACT 

We present carbon abundances of red giants in Milky Way globular clusters and dwarf spheroidal 
galaxies (dSphs). Our sample includes measurements of carbon abundances for 154 giants in the 
clusters NGC 2419, M68, and M15 and 398 giants in the dSphs Sculptor, Fornax, Ursa Minor, and 
Draco. This sample doubles the number of dSph stars with measurements of [C/Fe]. The [C/Fe] ratio 
in the clusters decreases with increasing luminosity above log (L/Lq) ~ 1.6, which can be explained 
by deep mixing in evolved giants. The same decrease is observed in dSphs, but the initial [C/Fe] of 
the dSph giants is not uniform. Stars in dSphs at lower metallicities have larger [C/Fe] ratios. We 
hypothesize that [C/Fe] (corrected to the initial carbon abundance) declines with increasing [Fe/H] 
due to the metallicity dependence of the carbon yield of asymptotic giant branch stars and due to 
the increasing importance of Type la supernovae at higher metallicities. We also identified 11 very 
carbon-rich giants (8 previously known) in three dSphs. However, our selection biases preclude a 
detailed comparison to the carbon-enhanced fraction of the Milky Way stellar halo. Nonetheless, the 
stars with [C/Fe] < +1 in dSphs follow a different [C/Fe] track with [Fe/H] than the halo stars. 
Specifically, [C/Fe] in dSphs begins to decline at lower [Fe/H] than in the halo. The difference in the 
metallicity of the [C/Fe] “knee” adds to the evidence from [a/Fe] distributions that the progenitors 
of the halo had a shorter timescale for chemical enrichment than the surviving dSphs. 

Subject headings: galaxies: dwarf- Local Group — galaxies: abundances — stars: evolution--stars: 
abundances 


1. INTRODUCTION 

Variations of elemental abundances within a stellar sys¬ 
tem are evidence for physical processes both internal and 
external to individual stars. Chemical evolution on a 
galactic scale leads to a dispersion in metallicity. As 
progressive generations of supernovae enrich the inter¬ 
stellar medium of a galaxy, stars that form later have 
higher metallicity than their older predecessors dTinslevl 
Il980l) . On the other hand, variations in certain light ele¬ 
ments can arise from internal stellar evolution. Gravita - 
tional settling and radiative levitation dBehr et al.lll999f ) 
deplete and enhance, respectively, the abundance of an 
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element in the photosphere of a star relative to the rest 
of the star. Nuclear burning at the base of a convec¬ 
tive zone that reaches the photosphere can also decrease 
or incr ease the observab le abundances of light elements, 
like Li (|Lind et all 1201 ll), C, N, O, and even s-process el¬ 
emen ts, like Ba and Rb (jLiang et al.ll200Ql: iKarakas et al.l 
I20H . ' This process is called astration when nuclear 
burning leads to the destruction of an element. 

Carbon exhibits abundance variations within a stellar 
population due to astration and galactic chemical evolu¬ 
tion. As a star ascends the red giant branch (RGB), 
it encounters episodes of mixing that bring material 
processed through the CNO cycle to the photosphere. 
Within 2-3 magnitudes of the tip of the RGB, C and 
O are depleted and N is enhanced relative to stars on 
the main sequence or the lower RGB. Proposed mech¬ 
anisms for mixing deep CNO-processed material into 
the upper atmosphere i nclude Rayleigh-Taylor instabil¬ 
ity llEggleton et al. l 20061) or , more likely, thermohaline 
mixing (jcharbonnel fe Zahnll2007l) . 

Globular clusters (GCs) are excellent sites to study the 
evolution of carbon. It is relatively easy to test whether 
variations in carbon abundances within a cluster are pri¬ 
mordial or modified by stellar evolutionary processes dur¬ 
ing the lifetime of the cluster because GCs are essentially 
single-age, mono-metallic populations at a uniform dis¬ 
tance. It is not necessary to control for age or metallicity. 
Carbon should not be astrated in metal-poor stars on the 
main sequence and subgiant branch, which do not have 
convective envelopes deep enough to dredge up CNO- 
processed material. Ho wever, GCs d o sho w primordial 
abund ance variations. ICohenl (|1999f ) and IBrilev et all 
(|2002f) found variations in [C/Fe] on the main sequences 
of M71 and M13, respectively. Because the photospheres 
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of main sequence stars are not affected by convection 
that can reach CNO-processed material, these carbon 
variations must be primordial. Regardless, astration on 
the upper RGB causes the most severe variations (e.g. , 
Sn^zefiyi^QSiySmith et^_aL| 119961 : iSmith fc Brilevl 120051 
200fiD . iMartell et alJ ( 2008cD showed that the rate of de¬ 
pletion depends_on meta lli city, as p redicted by thermo¬ 
haline mixing dCharbonnel fe Zahnll2007t) . 

Stars can also beco me very carbon-ri ch. “Ca rbon 
stars” (see the review by I Waller stein fc KnapplflOOSh dis¬ 
play extremely strong C 2 , CH, and/or CN bands, and 
they can be sig nifican tly redd er than the RGB in op- 
tical colors (e.g.. iBond fc Nef3119691 : iMould fe Aaronsonl 
Il979h . Carbon can form in helium-burning asymptotic 
giant branch (AGB) stars through the triple-a process. 
It can be brough t to the sur face of an AGB star in the 
third dredge-up fi benl 1975 1. This process is most ef¬ 
ficien t in stars around 2-2.5 M 0 (iMouhcine fc Lanconl 
120031) . which explains why most ancient GCs—where 
the hig hest-mass st ars are less than 1 M 0 —lack C stars 
dCohen et al.lll978l). C st ars also require at least 0.3 Gyr 
to form (IMouhcine 12002 ). which makes them good indi¬ 
cators of intermediate-age populations. 

Dwarf and RGB stars can acquire carbon through bi¬ 
nary mass transfer with a C star (see iLucatello et all 
120051: ICohen et al.ll2006l : IStarkenburg et al.ll2014D . C-rich 

RGB stars are rarely seen to have a photometric compan¬ 
ion because the AGB donor has almost always evolved 
into a faint white dwarf long ago. Instead, binarity of 
C-rich RGB stars is observed s pectrosco pically through 
varia bility in radial velocity ([McClure fe Woodsworthi 
I1990D . The recipient of carbon through mass transfer 
can exist in any evo lutionary stage, such as the main 
s equence (iDahn e t alJll977h. 

lAaronson fc Mouldl (|19 80) and lAaronson et all (119821 ) 
used the existence of carbon stars in dwarf spheroidal 
galaxies (dSphs) to argue for the presence of 
intermediate-age populations. These stars were selected 
for spectroscopy from their infrared colors, and their 
spectra showed very strong C 2 bands. However, car¬ 
bon enhancements need not be so large that th ey af¬ 
fect broadband colors. For example, iNorris et al. ( 201 0) 
discovered a C-rich ([C/Fe] = 2.3) red giant in the ultra¬ 
faint galaxy Segue 1. This star (Segue 1 7) could only 
have been found spectroscopically because its broadband 
colors are consistent with carbon-normal RGB stars. 

The statistics of carbon abundances in dwarf galaxies 
are less comp lete than foryGCs or t he Milky Way (MW). 
Although lAaronson fc Mouldl (119801) and their contem¬ 
poraries discovered carbon stars in dwarf gala xies, t hey 
did not quantify the carbon enhancement. iKinman et alJ 
( 1981 ) published in conference proceedings the first car¬ 
bon abundances from lo w-resol ution spectroscopy in a 
dSph, specifically Draco. iSmithl (119841 ) followed up with 
CH and CN indices in Draco, and iBelll (119851) pub¬ 
lished more low- resolution [C /H I abund ances in Draco 
and Ursa Minor. iShetrone et al.l (ll998bl ) later measured 
low-resolution [C/Fe] abundances for tw o stars in Sculp¬ 
tor, fin ding one to be mildly carbon-rich. iFulbrigh t et all 
( 2004 ) measured the first carbon abundance from high- 
resolution spectroscopy in a dSph, specifically Draco. 
lAbia et all (|2008l ) measured the 12 C/ 13 C isotopic ratio 
and the C/ O ratio—but not [C/H]— in t wo carbon stars 
in Carina. iColien fe Huand (120091 I2010D measured car¬ 


bon abundances in 18 stars in Draco and Ursa Minor 
from the G band of CH, and ISkuladottir et al.l (|2015l) 
measured C abundances for about 80 Sculptor giants 
from a near-infrared CN band. Recently, ultra-faint 
galaxies have received a lot of attention, includin g car¬ 
bon a bundance measurem ents in Segue 1 (]Norris et all 
l2010t IFrebel et al.l I2014D. Ur sa Major II and Coma 
Beren i ces (IFrebel et al.l l2010bl) . Bootes I (INorris et al. 
2010t lLai et al.l I20TTT and Bootes II dKoch fe Richl 


2014). Interestingly, all of these measurements are con¬ 


sistent with GC abun dances on the upper RG B except 
one star in Sculptor (ISkuladottir et al.l 12015 ) and four 
stars in Segue 1, including its three most metal-poor 
stars (INorris et al.l l2010h . The tendency of carbon-rich 
stars to appear amo ng the metal-poor population led 
IShetrone et, al.l (|2013[ ) to measure carbon ab undances for 
35 star s in the metal-poor dSph Dra co and IFrebel et all 
(I2010al) and IStarkenburg et all (120131 ) to explicitly target 
the most metal-poor stars in Sculptor. Still, Segue 1 and 
Sculptor contain the only examples of dwarf galaxy stars 
with quantified high carbon abundances. (As discussed 
above, C stars still exist in dSphs.) 

The MW halo contains a large number of C-rich stars, 
and their frequency increases with decreasing metallic¬ 
ity. However, the exact fraction of carbon-enhanced, 
metal-poor (C EMP) stars as a function of metallicity 
is contentious (Cohen et alJl2005t iMarsteller et alJl2005l: 


Lucatello et al. 120061: Frebel et all 120051: iCarollo et^l 

2014 lYong et, al.l 120131: iLee et al.) 120131 : IPI aero et al 


2014 ), partly because the search for low- met allicit y star s 


is fraught with selection biases (e.g JSchorck et all2009D . 
Very roughly, the fraction of stars with [Fe/H] < —2.5 
that are carbon-enhanced is about 20%. 

The MW halo is expected to be at least partly com¬ 
posed o f t he di srupted dSphs from past accretion events 
(jSearle fc Zinnl fl978t ). One apparent challenge to this 
idea was the discrepancy between th e [q/Fe] rat io of halo 
and dSph stars at [Fe/H] > —2 (IShetrone et al.lfl 998a. 
I2001L 20031 : IVenn et al.j l2004). However. iRobertson et al l 
(12005 ) and iFont et, all (I2006D showed that the discrep¬ 
ancy is in fact an expected outcome of hierarchical as¬ 
sembly. The inner halo, where most of the observed 
halo stars he, is composed of more massive satellites 
that do not resemble the present-day, surviving dSphs. 
At low enough metallicities, most elemental abundance 
ratios between dwarf galaxies and the MW halo agree 
(e.g., IFrebel et all 12010b ). though the neutron-capture 
elements are deficient in ultra-faint dwarf galaxies com- 
pared to th e halo, even at very low metallicities (e.g., 
IKoch et al.l[20l3l) . 

The fact that so few carbon-rich stars have been found 
in dSphs, even at low metallicities, might be troublesome 
for the theory of the hierarchical assembly of the MW 
halo. The fraction of CEMP stars in D raco , Ursa Mi¬ 
nor, and Sculptor is around 10% or less dCohen fc Huand 


2009tl20llilStarkenburg et all 20131 : IShetrone et al .1120131 


Skuladottir et al. nmlT However, this number does not 

include C-rich stars for which carbon abundances were 
not quantified. 

In this work, we add to the body of carbon abundance 
measurements in dSphs. We measured carbon abun¬ 
dances of 398 red giants (8 previously known) in the 
dSphs Sculptor, Fornax, Ursa Minor, and Draco. We also 
identified 11 giants that are too carbon-rich to quantify 
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their abundances. The s tars observed here are a subset 
of the stars observed by IKirbv et al.1 (|2009L 120 1 (1) , who 
measured Mg, Si, Ca, Ti, and Fe abundances in about 
3000 red giants in eight dSphs. 

2. SPECTROSCOPIC OBSERVATIONS 

We obtained Keck/DEIMOS ( Faber et al „ 2003 ) spec¬ 
tra of the carbon G band in red giants in MW GCs and 
dSphs. The GCs are NGC 2419, NGC 4590 (M68), and 
NGC 7078 (M15). The dSphs are Sculptor, Fornax, Ursa 
Minor, and Draco . We u sed the same dSph slitmasks as 
IKirbv et al.1 (|2009l [20Toh . who obtained redder DEIMOS 
spectra for measuring Mg, Si, Ca, Ti, and Fe abundances. 
We re-observed 3 of 5 slitmasks (198 of 376 member stars) 
in Sculptor, 1 of 5 slitmasks (125 of 675 members) in For¬ 
nax, 3 of 4 slitmasks (87 of 212 members) in Ursa Minor, 
and 3 of 8 slitmasks (181 of 298 members) in Draco. 

2.1. Target selection 

We selected stars for the DEIMOS slitmasks from 
a variety of photometric catalogs, which provided as¬ 
trometry, magnitudes, and co lors of the red giants that 
we observed. We used IStetsonl’s (12000H catalog for all 
three GCs. We supple mented IStetsonl’s list of M 68 stars 
with p hotometry from lAlcaino et al.l (|1990fl and iWalkeil 
(1199411 . We used Johnson V and Cousins I magnitudes 
for all GC stars. We corrected the magnitudes for red¬ 
dening and extinction based on E{B — V) values of 0.08, 
0.0 5, an d 0.10 for NGC 2419 M68, and M15, respectively 
(lllarridl 199(1 up date d 201(1^. 

IKirbv et al.l (I2009H desc r ibed the target selection for 
Sculptor, and IKirbv et al.1 ([2010) detailed the target se¬ 
lection for Fornax, Ursa M in or, an d Draco. For Sculp¬ 
tor, we used iWestfall et all ’s (1200611 photometric catalog 
in Washington M a n d Ti magnitudes. We used sources 
from iStetson et all s (119981 1 catalog in Cousins B and 
Joh nson V magnitud es for F ornax. For Ursa Minor, we 
used lBellazzini et~aIT s d2002h VI catalog, and for Draco, 
we used iSeeall et all ’s (I2007H gi catalog in the CFHT 
MegaCam system. 

We placed stars with colors and magnitudes appropri¬ 
ate for red giants on the slitmasks. For the GCs, we 
drew a polygon around the RGB in the color-magnitude 
diagram (CMD). Stars within the polygon were in¬ 
cluded in the slitmask. Brighter stars were preferred 
over fainter stars in the cases where slitmask design con¬ 
straints forced a choice among multiple red giants. For 
the dSphs, we selected stars from th e CM !) with the 
guidance of theoretical isochrones (see IKirbv et al.ll2010l 
for additional details). The Sculptor catalog addition¬ 
ally has DD021 magnitudes, which we used to discrimi¬ 
nate b etween re d giants and foreground dwarf stars (see 
IKirbv et al.ll2009j ) . 

2.2. Observations 

Table [T] lists the observation dates and exposure times 
of each slitmask. We used the 900ZD grating with a rul¬ 
ing of 900 lines mm -1 and a blaze wavelength of 5500 A. 
The grating was tilted such that the blaze wavelength 
fell on the center of the CCD mosaic. Slits were 0.7" 



rest wavelength (A) 

Figure 1. Portions of DEIMOS spectra around the G band for 
one star in each of the GCs and dSphs observed. The black line is 
the observed spectrum, and the red line is the best-fitting synthetic 
spectrum. The fit was performed in the green shaded region. The 
pale blue lines show synthetic spectra with [C/Fe] =b 0.2 different 
from the best-fitting value. 


wide. This configuration gives an approximate wave¬ 
length range of 4000-7400 A at a resolution of 2.1 A 
FWHM (i? ~ 2100 at 4300 A). The exact wavelength 
range of each spectrum depends on the location of the 
slit on the slitmask. The starting and ending wavelengths 
of the spectra vary by ~ 500 A across the slitmask. 

We obtained internal flat field and arc lamp exposures. 
We exposed a quartz lamp three times for 12 s and six 
times for 45 s, which is appropriate for the red and blue 
halves of the CCD mosaic, respectively. We also obtained 
two sets of arc lamp exposures. First, we obtained simul¬ 
taneous spectra of Ne, Ar, Kr, and Xe lamps for 1 s. This 
exposure is suitable to calibrate the wavelength solution 
for the red halves of the spectra. Then, we obtained 
a single exposure for several lamps. We controlled the 
amount of flux from each lamp by turning off lamps in 
the middle of the exposure. The lamps were Hg (1 s), Cd 
(5 s), Kr (6 s), Ar (12 s), and Zn (12 s). This exposure 
is appropriate for the blue halves of the spectra. 

All slitmasks were also observed with the 1200G grat¬ 
ing, which has a groove spacing of 1200 mm -1 and a 
blaze wavelength of 7760 A. The grating was centered at 
7800 A. Table 2 of IKirbv et al.l (120101) lists all of these ob¬ 
servations except for the n2419b and n4590a slitmasks. 
We observed n2419b with the 1200G grating and this 
7800 A setting on 2012 March 19 for 1800 s, and we 
observed n4590a on 2011 June 2 for 2400 s. 


ii 


http://www.physics.monaster.ca/$\sim$harris/mwgc.dat 


2.3. Reductions 
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Table 1 

DEIMOS Observations 


System 

Slitmask 

# targets 

Date 

Airmass 

Seeing 

(") 

Individual Exposures 

(s) 

Total Exposure Time 
00 

NGC 2419 

n2419b 

112 

2012 Mar 19 

i.i 

0.8 

3 x 900 

2700 

M68 

n4590a 

96 

2014 Feb 2 

1.6 

0.8 

1200 + 937 

2137 

M15 

7078d 

164 

2011 Jul 29 

1.0 

1.1 

3 x 600 

1800 


7078e 

167 

2011 Jul 30 

1.0 

0.9 

3 x 900 

2700 

Sculptor 

bscll 

86 

2011 Jul 31 

1.7 

0.8 

4 x 1200 + 900 

5700 


bscl2 

106 

2011 Aug 6 

1.8 

0.7 

2 x 1200 + 2 x 840 

4080 


bscl6 

91 

2011 Aug 4 

1.7 

0.8 

3 x 1200 + 1260 

4860 

Fornax 

bfor6 

169 

2011 Aug 5 

1.8 

1.3 

2 x 780 

6360 




2011 Aug 6 

1.8 

0.8 

1200 





2011 Aug 7 

1.9 

0.8 

3 x 1200 


Ursa Minor 

bumil 

125 

2011 Jul 29 

1.5 

0.7 

4 x 1200 + 600 

5400 


bumi2 

134 

2011 Jul 31 

1.7 

0.8 

4 x 1200 

4800 


bumi3 

137 

2011 Aug 4 

1.8 

0.6 

4 x 1200 

4800 

Draco 

bdral 

135 

2011 Jul 30 

1.4 

1.2 

5 x 1200 

6000 


bdra2 

144 

2011 Aug 7 

1.3 

0.7 

4 x 1200 

4800 


bdra3 

129 

2011 Aug 5 

1.3 

1.0 

5 x 1200 

6000 


We reduced the r aw data with the spec2d pipelin<3 
(iCooner e t al.1 120 121) devel oped by the DEEP2 survey 
team ( Newman et al.l 120131 ). This software is optimized 
for extended galaxies observed with the 1200G grating. 
We made moderate revisions to the pipeline to make it 
more suitable for bluer DEIMOS spectra. First, we mod¬ 
ified the software to accept separate flat lamp exposures 
for the red and blue CCDs. Second, we revised the arc 
lamp line lists for Hg, Cd, Kr, Ar, and Zn between 4158 A 
and 5496 A. We eliminated lines that we could not see 
in individual DEIMOS arc lamp exposures. Third, we 
changed the tolerances for automatic identification of 
arc lines so that the pipeline would include some of the 
weaker blue lines. Fourth, we changed the way that the 
two-dimensional spectrum is rectified. The curvature of 
the focal plane and the slit position angle can introduce a 
tilt in the spatial direction such that a single wavelength 
corresponds to a diagonal line rather than a single row of 
pixels. Previously, this tilt was removed by shifting CCD 
columns by integer values. We changed the code to allow 
non-integer shifts, which slightly increased the effective 
spectral resolution and signal-to-noise ratio (S/N). Fifth, 
we implemented a correction for atmospheric dispersion 
along the slit, which is important for our bluer spectra. 


Because light at 4000 A can be displaced along the slit by 
a few pixels from light at 7400 A, we curved the extrac¬ 
tion window to match the shape of the atmospheric dis¬ 
persion. Finally, we adopted the optimizati ons for point 
sources developed bv lSimon fc Ceha ( 200 7b 

The results of the reductions were wavelength- 
calibrated, sky-subtracted, one-dimensional spectra. 
Cosmic rays were removed before the one-dimensional 
spectra were extracted. The spec2d code also calculated 
and saved variance (error) spectra, which we used to es¬ 
timate random measurement uncertainty on [C/Fel (Sec¬ 
tion 13.4.1 ). 

Figure Q] shows one example spectrum from each stellar 
system in Table [Tj The figure shows spectra at a variety 
of S/N. Each panel lists the S/N, which is calculated as 
the median absolute deviation of the observed spectrum 
(black) from the best-fitting synthetic spectrum (red). 


12 http://www2.keck.hawaii.edu/inst/deimos/pipeline.html 
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Figure 2. The top panel shows a moderately carbon-rich AGB 
star in M15. The CN bandhead is at ~ 7850 A, and CN absorption 
is visible redward of the bandhead. The next two panels show CN 
absorption in two of the 11 carbon stars in the DEIMOS sample. 
The bottom panel shows a carbon-normal RGB star for compari¬ 
son. 


2.4. Carbon Stars 

We found 11 examples of carbon stars in our sample. 
The red spectra of these stars are dominated by very 
strong CN absorption, which precluded measurement of 
atmospheric parameters like [Fe/H]. The G bands are 
also strong enough that we would not have found our 
measurements of carbon abundances reliable due to sat¬ 
uration. Figure [2] shows examples of the near-infrared 
spectra of some of the se stars. We discuss carbon stars 
further in Sections 15.11 and [6] 
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3. SPECTROSCOPIC MEASUREMENTS 

Before we could measure carbon abundances from the 
spectra, we needed to measure the stars’ radial velocities 
(■ v r ) and atmospheric parameters. 

3.1. Radial Velocity Measurements 

We previously measured v r for all of the spectra 
from our red DEIMOS observations dKirbv et al.fI2009L 
I2010D . However, the wavelength calibration for the blue 
DEIMOS spectra is not as precise as the red spectra 
because of the relative sparsity of blue arc lamp lines. 
Rather than apply the known v r from the red spectra 
to the blue spectra, we measured v r directly from the 
G band in the blue spectra. Although tv,biue is not as 
precise as tv, re d (for example, for measuring galactic ve¬ 
locity dispersions), iv,biue is more suited to shifting the 
G band into the rest frame. 

We measured it^biue by cross-correlating the observed 
spectrum with a synthetic spectrum. The synthetic spec¬ 
trum had the atmospheric parameters previously mea¬ 
sured from the red spectra (see Section 15/21) and [C/Fe] = 
0.0. The spectra were synt hesi zed according to the pro¬ 
cedure described in Section l3Hl We performed the cross¬ 
correlation in the range 4200-4400 A, and we set v r to be 
the velocity corresponding to the cross-correlation peak. 
For the remaining analysis, we shifted the spectrum into 
the rest frame using this measurement of v r . 

We also evaluat ed the mem bership of each star by 
its radial velocity. iKirbv et~ahl ( 2010 . their Section 2.6) 
gave the details of the membership criteria. In summary, 
we eliminated stars more than three standard deviations 
away from the mean radial velocity. 


3.2. Atmospheric Parameters 

The measurement of carbon abundances (Section 13.41) 
requires knowledge of effective temperatures (T e ff), sur¬ 
face gravities (log g), metallicities ( [Fe/H]), and alpha en¬ 
hancements ([a/Fe]). IKirbv et~ahl ( 2000LT2(TT(ih measured 
these parameters for most of the spectra presented here. 
They obtained preliminary estimates of T e g and log g by 
comparing the stars’ colors and magnitudes to theoret¬ 
ical isochrones. They refined T e g and measured [Fe/H] 
and [a/Fe] by minimizing \ 2 between the observed red 
DEIMOS sp ectra and a large grid of synthetic spectra 
(IKirbvll20lif) . 

We adopted without modification the previously deter¬ 
mined atmospheric parameters. However, the DEIMOS 
slitmasks n2419b and n4590a were not included in 
the previous catalogs of atmospheric parameters. We 
measured these parameters in the same manner as 


IKirbv et al.l (120101) . 


3.3. Grid of Synthetic Spectra 

We measured [C/Fe] by comparing observed spectra of 
the G band to a grid of synthetic spectra spanning 4100- 
4500 A. We comp uted t his grid in a simil ar fa shion to 
the way in which lKirbv. Guliathakurta. fe Snedenl (|2008l ) 
computed their synthetic spectral grid spanning 6300- 
9100 A in order to measure [Fe/H] and [a/Fe] from red 
DEIMOS spectra. First, we compiled a line list from 
publicly available catalogs of atomic and molecular data. 
We refined the line list by comparing it to the spectra 


Table 2 

Line List 


Wavelength (A) a 

Species 

EP (eV) 

log 9f 

4100.005 

Mn II 

8.129 

-1.184 

4100.013 

NH 

2.202 

-3.421 

4100.014 

NH 

2.387 

-4.433 

4100.021 

NH 

2.335 

-4.289 

4100.022 

13 ch 

0.642 

-4.345 

4100.048 

NH 

2.824 

-2.929 

4100.051 

NH 

2.335 

-4.321 

4100.055 

NH 

2.599 

-2.932 

4100.084 

NH 

2.202 

-4.982 

4100.088 

Cri 

4.535 

-1.241 


References. — Atomi c data come from 
VALD dPiskuno v et al. l 119951; IKupka et al.l 119991) 
and NISTjEZamida ct al. 2014). CH data come 
from SCAN ( Jorgensen ct al. 1996). NH and CN 
data come from lKuruc A (I1992IT 

Note. — (This table is available in its entirety 
in a machine-readable form in the online journal. 
A portion is shown here for guidance regarding its 
form and content.) 
a Wavelength in air. 


of the Sun and Arcturus. Then, we computed MOOG 
synthetic spectra using ATLAS9 model atmospheres. 


3.3.1. Line List 


The G band has a high density of CH molecular ab¬ 
sorption lines. At the low resolution of our DEIMOS 
spectra, individual CH lines are indistinguishable from 
individual atomic and other molecular lines. Therefore, 
we needed to model the absorption due to many species, 
not just CH. We compiled a list of atomic and molecular 
transitions in the range 4100-4500 A. 

We started by including all the neutral and singly ion¬ 
ized atomic transitio ns listed in the Vienna Atomic Line 
Data basePl (VALD, iPiskunov et al.l 119951 : IKupka et al.1 
1999) in our spectral range of interest. We discarded 
transitions with excitation potentials (EP) greater than 
10 eV or with oscillator strengths log gf < —5. 

We then compiled a separate atomic line list from the 
National Institute of Standards and Technology (NIST) 


databasco (jKramida et al.l 120141) . This database is 
very incomplete, but it generally has more trustworthy, 
laboratory-measured atomic transition probabilities than 
other databases. We replaced entries in the VALD line 
list with NIST entries where available. 

Next, we added molecula r transit io ns. W e used the 
CN and NH line lists of IKurucd ( 1992TP’! and the 
SCAN line lislP^l for CH from iJoreensen et al.l (I1996D . 
(|Masseron et al.l 120 141 published a new CH line list long 
after we computed the spectral grid and shortly before 
we submitted this article.) We preserved the isotopic in¬ 
formation for each CH transition so that we could change 
12 C/ 13 C in the spectral syntheses. 

No line list is perfectly accurate or perfectly com¬ 
plete. We tested our own line list against observed, high- 


13 http://vald.astro.univie.ac.at/$\sim$vald/php/vald.php 

14 http://www.nist.gov/pml/data/asd.cfm 
lo http://kurucz.harvard.edu/ 

16 http:// tow. astro.ku.dk/$\sim$uffegj/ 
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Table 3 

Grid of Synthetic Spectra 


Parameter 

Start 

End 

Step 

A 

4100 A 

4500 A 

0.02 A 

T e fi 

3500 K 

5600 K 

100 K 


5600 K 

6400 K 

200 K 

log 9 

0.0 

4.0 

0.5 

[Fe/H] 

-4.0 

0.0 

0.1 

[cx/Fe] 

-0.8 

+ 1.2 

0.1 

[C/Fe] 

-2.4 

+ 1.0 

0.2 


+1.0 

+3.0 

0.4 


+3.0 

+3.5 

0.5 


resolu tion spectra of the Sun and Arcturus (iHinkle et al.l 
2000). We co mputed synthe tic spectra for both stars us¬ 
ing MOOCH (|Snedenl|1973f) with an und ated treatment 
of electron scattering (iSobeck et al. I201lh . For the Sun, 
we used an ATLAS9 (|Kurucdll993l) model atmosphere 
with T R ff = 5798 K a nd log t? = 4.44. For Arcturus, we 
used iPeterson et ahl ’s ()1993l ) model atmosphere. Where 
atomic lines did not match the observed spectra, we 
changed oscillator strengths. We also introduced a hand¬ 
ful of artificial Fe I lines to mimic observed lines that did 
not appear in our original line list. We avoided chang¬ 
ing molecular line strengths. Table [5] gives the final line 
list that we us ed in our computation of synthetic spectra 
(Section 13.3.21) . 

3.3.2. Spectral Synthesis 

We used MOOG with ou r ex isting grid of ATLAS9 
model atmospheres dKirbvl [201 lh to compute synthetic 
spectra at a range of atmospheric parameters. Table [3] 
shows the parameters of the spectral grid: the lower 
and upper ranges and step sizes of wavelength, effective 
temperature, surface gravity, metallicity, alpha enhance¬ 
ment, and carbon abundance. We used 256 CPUs to 
compute the 4,835,376 spectra. Although [C/Fe] varied 
in the spectral syntheses, the input model atmospheres 
were all computed assuming [C/Fe] = 0. 

As stars ascend the RGB, they dredge up products 
of the CNO cycle. We examine this phenomenon in 
Section 01 but many others h ave investigated it be¬ 
fore. In particular, iKeller et al.l (l2001f ) showed how the 
12 C/ 13 C isotope ratio varies as a function of evolution¬ 
ary state. The following relation, which we deduced 
from lKeller et all ’s Figure 4, approximates how the ratio 
changes with decreasing surface gravity: 

12 C/ 13 C = 50 if logg > 2.7 

12 C/ 13 C = 63 log<? — 120 if 2.0 < log g < 2.7 m 

12 C/ 13 C = 6 if log< 2.0 (1) 

We used Equation [T] in computing the grid of synthetic 
spectra. 

3.4. Carbon Abundance Measurements 

We measured carbon abundances by comparing the ob¬ 
served spectra to our new grid. This process involved 
several steps. 

1 ' http://www.as.utexas.edu/$\sim$chris/moog.html 


1. Continuum division: MOOG generates synthetic 
spectra normalized to the continuum. We divided 
each observed spectrum by a continuum fit so that 
it could be compared to the synthetic spectra. We 
fit the spectrum with a spline with a breakpoint 
spacing of 200 A. Pixels with fluxes that were 
greater than 5 standard deviations above or 0.1 
standard deviations below the spline were itera¬ 
tively excluded from the fit so that stellar absorp¬ 
tion was not included in the continuum determi¬ 
nation. The spectrum was divided by the spline. 
This continuum was refined in step 3. 

2. Preliminary [C/Fe] measurement: We searched the 
spectral grid for the best-fitting value of [C/Fe]. 
The synthetic spectra were smoothed through a 
Gaussian kernel (2.1 A FWHM) to match the ob¬ 
served resolution. We fixed T e g, logg, [Fe/H], and 
[a/Fe] at the values previously determined from the 
red DEIMOS spectra (see Section 13131) . We used 
the Leven berg-Marquardt optimization code MP- 
FIT (!Markwardtll2012lf to find the value of [C/Fe] 
that minimized y 2 in the spectral range 4260- 
4325 A (green region in Figures [1] and [3]), where 
the spectrum is most sensitive to changes in car¬ 
bon abundance. 

3. Continuum refinement: The continuum determina¬ 
tion in step 1 does not take into account the known 
variation in flux as a function of wavelength. We re¬ 
fined the continuum by dividing the observed spec¬ 
trum by the best-fitting synthetic spectrum deter¬ 
mined in step 2. We fit that quotient with a spline 
with a breakpoint spacing of 150 A. This spacing 
is slightly finer than in step 1 for a slightly more 
detailed continuum determination. We iteratively 
discarded pixels that were above or below the fit 
by more than one standard deviation. 

4. Refined [C/Fe] measurement: We repeated steps 2 
and 3 until [C/Fe] changed by less than 0.001 be¬ 
tween iterations. 

In some cases, the spectra allowed only measurements 
of upper limits. We computed y 2 contours for each star 
at several values of [C/Fe] around the minimum y 2 . For 
stars whose y 2 contours did not rise by at least 1 on 
both sides of the minimum, we calculated a 2cr upper 
limit. The upper limit was the value of [C/Fe] at which 
y 2 was 4 above the minimum y 2 . 

FigureQ]shows model spectra of the best-fitting G band 
spectrum in red. Also shown in light blue are model spec¬ 
tra deviant from the best fit at [C/Fe] ± 0.2. Figure [3] 
shows two pairs of spectra of stars—one pair in Sculp¬ 
tor and one pair in Ursa Minor— with similar luminosity 
and metallicity but different carbon abundances. The fig¬ 
ure demonstrates that the stars measured to have larger 
[C/Fe] have noticeably stronger G bands. 

We also computed the S 2 (CH) index (iMartell et all 
l2008a|) . This index is defined to have minimal sensitiv¬ 
ity to nitrogen abundance. It is the ratio of flux in the 
G band to flux in two side bands. The blue side band 
starts at 4212 A. We measured S^CH) only in spec¬ 
tra with minimum rest wavelengths less than 4212 A. 
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4260 4280 4300 _ 4320 4340 

rest wavelength (A) 


Figure 3. G band spectra of pairs of stars with similar luminos¬ 
ity, similar metallicity, and differing carbon abundance in Sculptor 
(top) and Ursa Minor (bottom). The stellar parameters listed on 
the left (right) of each panel correspond to the solid black (dashed 
blue) spectrum. As in Figure [T] only the green shaded region is 
used in the measurement of [C/Fe]. 

iMartell et al.1 (l2008af ) defined the band to be measured 
in units of ADU per pixel. However, the throughput of 
DEIMOS in the vicinity of the G band diminishes rapidly 
toward the blue. Consequently, we measured S 2 (CH) in 
continuum-corrected flux rather than the original ADU 
per pixel. 

Table S] gives the [C/Fe] and S^CH) measurements 
for our GC and dSph sample. The last column iden¬ 
tifies the 11 carbon stars for which we were unable 
to measure [Fe/H] or [C/Fe]. The table also gives 
photometric information, including the absolute lumi¬ 
nosity for each star relative to the solar luminosity. 
We used b olometric correc t ions b ased on theoretica l 
isochrones (jDemaraue et al.1 12004 iGirardi et al.l I2004D 
and previously measured distance moduli t o co mpute 
the luminosities: 19.83 for NGC 2419 (lHarris et, al.1 
Fl99l . 15 21 for M68 (IMcClure et al.l 11 98711 . 15.08 
for M15 (iDurrell fc Harris! I1993H . 19.67 for Sculptor 
(iPietrzvriski et al.l I2008D , 2 0.72 for Fornax ( Rizz i et al.l 
1200711 . 19.18 for Ursa Minor (IMighell fc Burkd19990 . and 
19.84 for Draco dBellazzini et al.l 1200211 . The table also 
includes T e g, log g, [a/Fe], and S/N, calculated as the 
median absolute deviation of the observed spectrum from 
the best-fitting spectrum in the range 4250-4350 A. 

3.4.1. Estimation of Uncertainties 

Although MPFIT provides a covariance matrix for 
each minimization, this estimate of uncertainty on [C/Fe] 
does not account for uncertainty in the placement of the 
continuum. Instead of using the covariance matrix, we 
calculated Monte Carlo uncertainties on [C/Fe] by re¬ 
sampling the observed spectrum. 

After we determined the best-fitting carbon abun¬ 
dance, we added noise to the spectrum after it was 
continuum-normalized according to step 1 in Section [3.41 
The new flux in each noise-added pixel was the original 
flux plus a value drawn from a Gaussian random distri- 



Figure 4. The response of [C/Fe] to changing [a/Fe] in the stellar 
atmosphere. The a:-axis shows the difference between the [a/Fe] 
values used in our fiducial measurements of [C/Fe] and [a/Fe] = 
0 . 0 . 

bution with a standard deviation equal to spec2d’s esti¬ 
mate of the flux error in that pixel. We repeated steps 
2-4 above with the noise-added spectrum. 

We repeated this process for a total of 1000 noise re¬ 
alizations. We took the base uncertainty on [C/Fe] as 
the standard deviation of all 1000 trials. Although this 
procedure accounts for random uncertainty, it does not 
include all sources of systematic error. We added 0.1 dex 
as an estimate of systematic error in quadrature with 
the Monte Carlo uncertainty for the final uncertainty on 
[C/Fe]. This systematic uncertainty is approximately the 
value that is required to account for the rms of the dis¬ 
tribution of our measurements of [C/Fe] versus literature 
measurements (Figure [5] discussed in Section 13.4.31) . 

We computed errors on ^(CH) from the same noise- 
added spectra. We did not add any systematic error to 
these measurements. 

3.4.2. Effect of Nitrogen and Oxygen 

The measurement of carbon could be affected by the 
abundances of other elements. First, other elements that 
form diatoms with carbon affect the molecular equilib¬ 
rium, changing the amount of carbon available to form 
CH. Second, some elements contribute to the opacity in 
the stellar atmosphere. Changes in their abundances can 
affect the structure of the atmosphere and the strength 
of the G band. The most important elements to consider 
for these effects are nitrogen and oxygen because they are 
abundant, and they form diatoms with carbon. Unfortu¬ 
nately, we can measure neither element directly from our 
spectra. However, we can quantify their possible effects 
on our carbon abundance measurements. 

We explored the effects of molecular equilibrium be¬ 
tween CH and CN by making syntheses with MOOG 
with varying abundances of N. We took three stars in 
M15 ([Fe/H] ss —2.4) as test cases: 13701, which is be¬ 
low the RGB luminosity function bump; 43026, which 
is at the top of the RGB; and 47866, which is carbon- 
enhanced (see Section SI). We also tested two stars in 
Fornax to explore more metal-rich ([Fe/H] « —1) stars: 
52752 at log(L/Lg) = 2.0 and 40828 at the top of the 
RGB. We synthesized spectra with [N/Fe] values of —2.0, 
-1.0, -0.5, -0.2, 0.0, +0.2, +0.5, +1.0, and +2.0. In no 
case did the flux of any pixel in the G band change by 
more than 1%. 

Our carbon abundance measurements already account 
for varying elemental abundances that might affect the 
structure of the stellar atmospheres. We previously mea¬ 
sured four a elements (Mg, Si, Ca, and Ti) from red 
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DEIMOS spectra (iKirbv et al.ll2010l) . The stellar atmo¬ 
spheres we used to measure carbon abundanc es reflect 
those [a/Fe] ratios. IVaudenBerg et all (|2012h showed 
that at temperatures found in the atmospheres of cool 
giants, Mg and Si are the most important metals for de¬ 
termining the opacity. As a result, changing [a/Fe] (in¬ 
cluding [Mg/Fe] and [Si/Fe]) can alter the atmospheric 
structure and therefore affect the measurements of all 
elemental abundances. 

We recomputed all of the [C/Fe] measurements assum¬ 
ing [a/Fe] = 0.0 rather than the values previously deter¬ 
mined bv IKirbv et all ( 2010 ). Here, a refers to O, Ne, 
Mg, Si, S, Ar, Ca, and Ti. This value of [a/Fe] is used 
to compute both the model atmosphere and the spectral 
synthesis. Hence, the syntheses include the effects of the 
varying atmospheric structure and the changing molec¬ 
ular equilibrium between CH and CO. Figure 0] shows 
the response of [C/Fe] to the change in [a/Fe]. Changing 
[a/Fe] by ±0.6 dex—much more than the mean uncer¬ 
tainty in [a/Fe], 0.2 dex—results in a change of approx¬ 
imately ±0.1 dex in [C/Fe]. Forcing [a/Fe] to be 0.0 
results in |A[C/Fe]| <0.1 dex for 94% of the dSph stars 
and | A [C/Fe]] < 0.05 dex for 75% of the dSph stars. We 
conclude that uncertainty in [a/Fe] is not a major source 
of uncertainty in our carbon abundance measurements. 

In conclusion, uncertainties in nitrogen, oxygen, and 
other a element abundances are minor contributors to 
the error budget of our [C/Fe] measurements. 

3.4.3. Validation 

We compared our low-resolution measurements of 
[C/Fe] with previous high- and low-resolution measure¬ 
ments of [C/Fe] for the same stars. We found overlapping 
stars with high- resolution measureme nts in t he literature 
for N GC 2419 (|Cohgn_j{%Kjrbyt l2012fl. M68 dVenn et al. 
2012 ), Draco dFulbright et all 120041 ICohen fc Huang 
2009), and Ursa Minor (ICohen fc Huangl I2010D . 
Shetrone et al.l (120131) also published measurements of 
[C/Fe] in Draco based on low-resolution Keck/LRIS spec¬ 
troscopy. All of these measurements are also based on the 
G band. 

The left panel of Figure [5] shows the comparison 
between [C/H] measured from DEIMOS and [C/H] 
as published in the aforementioned references. The 
right panel shows [C/Fe]. We shifted all abundances 
to ou r adopted sol ar abundance scalcQ A(C) = 
8.56 (lAnders fc Grevessel I1989D and A(Fe) = 7.52 
(lSneden et al.lll99^ ). The deviations reach up to 0.8 dex. 
The deviations are also about the same for [C/Fe] as 
they are for [C/H], which indicates that differences in 
atmospheric parameters like T e g- are not the reason for 
the discrepancy. On the other hand, the degree of de¬ 
viation depends largely on the source of the measure¬ 
ment. For exa mple, our [ C/H] measurements are 0.3- 
0.4 dex above [Venn et all s [C/H] measurements from 
VLT/FLAMES, b ut our measurem en ts are on averag e 
~ 0 .4 dex below ICohen fc HuangT s. ICohen fc Kirbvf s. 
and iFulbright et all s measurements fro m Keck/H IRES. 
Our measurements largely agree with iShetrone et all ’s 
Keck/LRIS measurements. 

18 A(X) = 12 + log(n(X)/n(H)) where n(X) is the number den¬ 
sity of element X. 


The origin of the discrepancies might be spectral res¬ 
olution or the details of the carbon abundance measure¬ 
ments. At first glance, the spectral resolution is suspect 
because our low-resolution meas urements agree most 
with another low-resolution study ([Shetrone et_alj 20j_5 ). 
However, our measurements disagree with [Venn et all ’s 
(I2012D in the opposite direction from the other two high- 
resolution studies. Therefore, we investigated the possi¬ 
bility that differences in the line list caused the offset. We 
compared our line list (Table [2]) with J. Cohen’s line list. 
Between 4274 A and 4232 A, our list contains 2.5 times 
more CH lines. The higher density of lines makes the 
G band stronger. Hence, for a fixed observed G band 
strength, we would measure a lower carbon abundance 
with our line list than with Cohen’s line l ist. Our line 
list has a similar density of lines to IShetrone et all ’s list. 
Therefore, it is plausible that a different line list is the 
main source of disagreement of carbon abundance mea¬ 
surements from the G band. 

4. ASTRATION OF CARBON ON THE UPPER RGB 

The destruction of carbon on the RGB above the lumi¬ 
nosity function bump is well documented (e.g. , iSuntzeffl 
119811: ICarbon et ahlll982l : ISmith fc BrilevH2006l ). As pre¬ 
vious investigators have done, we searched for signatures 
of carbon astration in GCs. GCs are the easiest stel¬ 
lar systems to isolate stellar luminosity as an indepen¬ 
dent variable because the cluster stars have similar ages, 
metallicities, and distances. 

Figure [6] shows the trend of carbon enhancement 
([C/Fe]) with luminosity. The carbon abundance re¬ 
mains flat until giants reach log (L/Lq) = 1.6, where¬ 
upon [C/Fe] declines, presumably due to mixing that be¬ 
gins at the luminosity function bump. We fit a line to the 
combined trend of [C/Fe] with [Fe/H] for all three GCs, 
fixing the luminosity of the bump at log(L/L 0 ) = 1.6. 
The best fit line, represented in green in Figure |6l is 

[C/Fe] = 0.12 ± 0.06 if log(L/L 0 ) < 1.6 

[C/Fe] = (1.42 ± 0.06) ( . 

-(0.82 ± 0.02) log(L/L 0 ) if log(L/L 0 ) > 1.6^' 

The amount of carbon astration in GCs could be com¬ 
puted from this empirical trend, but we also considered 
t heoretical es timates of the amount of carbon depletion. 
iPlacco et al.l (I2014D calculated corrections based on mod¬ 
els of mixing on the upper RGB. The corrections de¬ 
pend on surface gravity, metallicity, and uncorrected car¬ 
bon abundance Fu We used their models to compute a 
correction for each star in our sample. Figure |6] shows 
in faded colors IPlacco et all ’s corrections applied to the 
GC sample. The faded points have less dependence on 
luminosity than the uncorrected points. Although the 
empirical li near fit to the data (Equation [2]) is a better 
fit than IPlacco et aP ’s models, the models are ab initio 

19 The corr ectio n depends slightly on [N/Fe]. We assumed 
[N/Fe] = 0. IGratton et alj (120001 ) showed that [N/Fe] climbs 
to about 0.5 at the tip of the RGB in mctal-poor field stars. 
IPlacco et all 's < 20141 1 corrections change by < 0.05 dex for changes 
in [N/Fe] of 0.5 dex. For the range of stellar parameters in our sam¬ 
ple, the difference from the [N/Fe] = 0 correction exceeds 0.1 dex 
only for [N/Fe] > +1, a condition that very few, if any, of the stars 
in our sample would satisfy. 
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Figure 5. Comparison between Keck/DEIMOS and literature carbon abundances for stars in the GCs NGC 2419 and M68 a nd the 
dSphs Draco and Urs a Mi nor. All o f the literature measurements come from Keck/HIRES except for the studies of iVenn et al.l 1 12013 . 
Magellan/MIKE) and lShetrone et al. l 1 2001] , Keck/LRIS). 


predictions. They can be applied to data in environ¬ 
ments that are not as tig htly controlled as GCs. We 
use the iPlacco et al.1 (|2014r) corrections for dSph stars in 
Section [5] 

IMartell et al.l (|2008cT ) found that the rate of carbon 
depletion as giants evolve depends on metallicity, which 
is consistent with the thermohaline mixing mechanism 
for the_change_ i n elemen ta l composit ion with luminos¬ 
ity ( Charbonnel fc Zahnl 120071 ICharbonnel fe Lagardel 
l2010f) . Unfortunately, we can not test for metallicity 
dependence because the GC s we ob s erved have similar 
metallicities. According to IHarrisI (U996|l. the metal- 
licitie s are [Fe/H] = —2.15 (NGC 2 419, TSuntzeff et al.1 
19881). -2.28(M68 . iLeeetaLl [Mil , and -2.37 (M15, 
Takeda et al.l 120091) . IMartell et al.1 found carbon deple¬ 
tion rates of d[C/Fe]/dlog(i/LQ) = —0.7 to —1.5 above 
the luminosity function bump for the metallicity range of 
these three GCs. Our value of —0.82 ±0.02 (Equation[2]) 
falls in this range. 

GCs have prim ordial abundance var iations 
dGratton et al.l 120041) . including carbon (e.g., iCohenl 
I1999D . The dispersion in [C/Fe] at fixed luminosity can 
give a sense of the differential enrichment in carbon 
between the first and subsequent generations of star 
formation. We calculated this dispersion by finding the 
residuals between [C/Fe] and the kinked linear trend. 
The trend is similar to Equation 0 but we fit each GC 


individually to calculate the dispersion. We then nor¬ 
malized the residuals by the measurement uncertainty: 
A = ([C/Fe] — [C/Fe](trend))/<5[C/Fe]. If there is no 
intrinsic scatter in [C/Fe] beyond the measurement 
uncertainty, then the standard deviation of A should be 
1. Instead, we found 1.6 for NGC 2419, 1.5 for M68, and 
2.1 for M15. (We excluded star 47866, discussed below, 
from M15.) This measurement of dispersion demands 
that we estimated measurement uncertainties accu¬ 
rately. If we underestimated the uncertainties, then the 
standard deviation of A will appear erroneously larger 
than 1. Taken at face value, we have found evidence 
for intrinsic scatter in [C/Fe]. We see no evidence for 
bimodality in [C/Fe], but the scatter is not much larger 
than the uncertainties. No netheless, CN bands d o show 
bimodality in some GCs dSmith fc Brilevl 120051 120061 : 
IMartell et al.1 l200~8bl) . reflecting that primordial C and 
N abundances in clusters were bimodal, independent 
of the present-day astration on the upper RGB. This 
bimodality can even be found on the main sequenc e of 
s ome clusters ( e. g., 47 T ucanae. lHarbeck et al.ll2003l) . 

iCohen et al.1 (201 111. ICohen fe Kirbvl ( 2012D . and 
iMucciarelli et ah ( 20121) discovered that NGC 2419 has 
an extremely unusual chemical composition. The Mg 
abundances vary from star to star by a factor of 25, 
and the variation is anticorrelated with abundances of 
other elements, including K and Ca. In a sample of 13 
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Figure 6. The decrease in carbon enhancement with red giant lu¬ 
minosity in two globular clusters. The giants in NGC 2419 are sep¬ 
arated into calcium-enhanced (red triangles) and normal-calcium 
groups (red circles). The green lines illustrate of the destruction 
of carbon as photospheric material is processed by the CNO cy¬ 
cle on the upper RGB (Equation [2]) • The faded points show the 
[C/Fe] measurements corrected for the destruction o f carbon with 
i ncrea sing luminosity according to the calculations of IPlacco et al.l 

m3)- 

red g iants with high-resolution spectra, I Cohen fe Kirbvl 
(2012:) did not find any correlation between C and any 
of the elements with an intracluster dispersion. We 
looked in our larger sample for such a correlation by 
separating the NGC 2419 stars into Ca-enhanced and 
Ca-normal. We estimated Ca enhancement by measur¬ 
ing the equivalent widths (EWs) of the near-infrared Ca 
triple t. We converted th es e EW s to metallicity estimates 
using lStarkenburg et ahl ’s (|2010l) formula. The stars with 
[Fe/H] CaT < —2.2 (Ca-normal) are represented by blue 
circles in Figure [6] Stars with [Fe/H] CaT > —2.2 (Ca- 
enhanced) are red triangles. There is no noticeable dif¬ 
ference between the two populations. Hence, we confirm 
that carbon did not participate in the unusual nucleosyn¬ 
thesis that affected heavier elements in NGC 2419. 

Figure [7] shows an alternate way to display the destruc¬ 
tion of carbon as red giants ascend the RGB. The CMD 
color-codes the stars, with cooler colors corresponding to 
lower [C/Fe]. In general, the plotting color is cooler at 
the top of the RGB compared to less luminous stars. In 
M68 and M15, which are the GCs observed to sufficient 
depth, most of the stars below the luminosity function 
bump ( Iq ~ 14) are colored green. This color corresponds 
to the constant pre-bump value of [C/Fe] = +0.12. 


Note the one carbon-rich star in M15, which is also 
visible in Figure [6] Star 47866, a confirmed radial ve¬ 
locity member of M15, stands out as lone orange point 
([C/Fe] = +0.88 ± 0.10) in Figure [7] in a sea of dark 
green points ([C/Fe] ~ —0.2). The carbon abundance is 
high enough that the near-infrared spectrum shows some 
CN absorption (Figure [7]). This star was also previously 
identified by J. Cohen to be carbon-rich from an unpub¬ 
lished Keck/HIRES spectrum. That spectrum shows a 
weak C 2 bandhead at 5163 A in addition to a very strong 
G band. Also note that this star lies slightly blueward of 
the main RGB locus. The CMD position indicates that 
star 47866 is an AGB star. 

Car bon- enhanced s t ars in GCs are rare. iCohen et al.l 
lfl97j and lMouhcinel (|2002t l explain that C stars appear 
after several hundred million years, long after the end of 
star formation in a GC. None of the GC stars would have 
been pre-enriched by carbon from an AGB star. Addi¬ 
tionally, the lowest-mass AGB star that produces copious 
amounts of carbon is about 2 M 0 . Such a star would live 
for about 1 Gyr. Hence, all carbon-producing AGB stars 
in GCs should have died long ago. Nonetheless, there are 
a few examples of C-rich stars in GCs. T he GC 1 0 Cen- 
tauri hosts at l east seven k n own C stars (Harding] 1962; 
Dickensl 119721: iBondl 119751 : iCowlev fe Cramptonl 1985 : 
van Loon et al.ll2007f) . However, uj Centauri is known to 
have a n unusually long star formation history (SFH) for 
a GC ([Hilker et aLll2004l iVillanova et al. 1120071 ). Regard¬ 
less, C st ars can be found in even m ore ordinary GCs , 
like M14 dCote et al.lIl997h . Lvng a 7 dMatsunagal 120061 : 
iFeast et al.ll2013fl . and NGC 6426 (|Sharina et al.ll2012D . 

C stars can be formed as the result of mass trans¬ 
fer from an AGB star. In most cases, the AGB star 
has now evo lved into a w hite dwarf. There is pho¬ 
tometric (iBohm-Vitense et al.l l2000fl and spe ctroscopic 
(jLucatello et al.ll2005b iStarkenburg et al.ll2014f) evidence 
that most or all C stars—especially those enhanced in s- 
process elements—have binary companions. Binarity is 
often diagnosed through radial velocity variability over 
multiple measurements. Although we observed 47866 
twice, the observations were separated by only 24 hours. 
We found no evidence for binarity from this very short 
time baseline. An alternative possibility is that the ex¬ 
isting C-rich AGB star formed in a stellar merger, the 
result of which would be mas sive enough t o gen erate car¬ 
bon in its AGB phase le.g.. IFeast et al.l 1201 31). Stellar 
mergers are a plau s ible origin for bl ue stragglers in GCs 
(|Mateo et al.lll990l : iSills et aT1l2013T) . so it is conceivable 
that this carbon-rich star is also the product of a merger. 

5. CARBON IN DWARF SPHEROIDAL GALAXIES 

Unlike GCs, dSphs contain stars with a range of ages 
and metallicities. The more complex stellar populations 
complicate the evolution of carbon. First, dSph stars 
need not start their lives with the same C abundance. 
GC stars were formed from well-mixed gas polluted only 
by Type II supernovae (SNe) and AGB stars in a limited 
mass range. Figure [6] and Equation [2] show that all of the 
pre-bunrp stars within a GC have roughly the same value 
of [C/Fe] ~ +0.12. On the other hand, the progenitors 
of dSph stars can have a range of metallicities. They can 
also form from inhomogeneously mixed gas polluted by 
SNe of various types and AGB stars of various masses 
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Figure 7. Color—magnitude diagrams for the GCs and dSphs. Point color corresponds to [C/Fe]. Hollow circles are spectroscopic targets 
that are non-members or spectra where [C/Fe] was unable to be measured. Red, five-pointed symbols indicate very carbon rich stars for 
which we did not attempt to measure [C/Fe]. 


and metallicities. Second, the star formation duration 
of dSphs is long enough for 2-5 Mq AGB stars—the 
ones that generate copious amount of carbon—to pollute 
the interstellar medium while the galaxy is still forming 
stars. As a result, AGB stars influence the evolution 
of carbon differently in dSphs than in GCs. Third, the 
stellar density of dSphs is significantly lower than GCs, 
which reduces the stellar merger rate. Consequently, the 
merger formation channel for prematurely C-rich stars 
(possibly like 47866 in M15) is suppressed. 

The bottom panels of Figure[7]show that the decline of 
[C/Fe] with increasing stellar luminosity is just as present 
in dSphs as in GCs. However, there is a larger spread in 
[C/Fe] at fixed luminosity in the dSphs. Figure [5] shows 
the astration of carbon in dSphs, just as Figure [5] shows 
the same for GCs. The green lines in Figure [5] are the 
same as in Figure [6j They were not re-fit to the dSph 
data. The slope of the line roughly tracks the decline of 


[C/Fe] in the dSphs, especially in Sculptor and Draco. 

The scatter in [C/Fe] at fixed luminosity exceeds the 
scatter expected from measurement uncertainty. The 
residuals are correlated with [Fe/H], which is represented 
by the color coding in Figure [8] The residual metallicity 
trend is especially apparent in Ursa Minor and Draco, 
the most metal-poor of the four dSphs we observed. The 
stars with the lowest [Fe/H] have the highest [C/Fe] at 
fixed luminosity. This trend indicates that C abundance 
did not increase as quickly as Fe abundance as star forma¬ 
tion progressed in the dSphs. Later forming stars, which 
have higher [Fe/H], started with lower [C/Fe] ratios than 
their predecessors. 

Figure [9] shows the trend in another way. The x-axis is 
[Fe/H] instead of luminosity. Instead of [C/Fe], we plot 
[C/Fe] corr , which is corrected for the depletion due to 
astrat ion on the upper RGB according to iPlacco et alJ 
(|2014l) . This correction removes the luminosity depen- 
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Figure 8. Carbon enhancement as a function of stellar luminosity in dSph stars. The green line ('Equation [2}, copied from the GC data 
in Figure [6] shows the destruction of carbon as photosph eric materia l i s proc essed through the CNO cycle on the upper RGB. The color 
coding corresponds to iron abundance. The dotted line is lAoki et all s (|20071 ) dividing line for carbon enhancement. Upward-pointing red 
arrows indicate very carbon rich stars for which we did not attempt to measure [C/Fe]. Downward-pointing gray arrows are 2 cr upper 
limits. 
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Figure 9. Carbon enhancement corrected for luminosity accordin g to | Placco et ahl (120141) as a function of stellar metallicity. The color 
coding corresponds to stellar luminosity. The dotted line is one of IPlacco et all s definitions for carbon enhancement. Downward-pointing 
gray arrows are 2cr upper limits. The C stars are not plotted because we did not measure [Fe/H] for those stars. 
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Figure 10. Carbon e n hance ment corrected for luminosity ac¬ 
cording to IPlacco et al.l (|2014f ) for all four dSphs. The theoret¬ 
ically predicted, IMF-averaged yields for nucleosy nthetic sources 
are also show n: AGB stars (solid line, [Karakas 2010J), Type II SNe 
(dashed line, Nomoto et al.ll2006h . and Type la SNe (dotted line, 
llwamoto et a .1119991 1. The AGB and Type II SN yields depend 
on meta llicity. AGB stars could be active at [Fe/H] < —2.3, but 
IKarakasl ’s models are not available at lower [Fe/H]. Six-pointed 
symbols indicate the metallicities at which the yields were com¬ 
puted. A decline in [C/Fe] corr can be explained by the metallicity 
dependence of AGB yields combined with a growing contribution 
of Type la SNe. 


dence so that we may isolate [Fe/H] as the indepen¬ 
dent variable. Figure [§] confirms that [C/Fe] corr declines 
with increasing [Fe/H] in Draco. The trend is also ap¬ 
parent in Ursa Minor and Sculptor, where stars with 
[Fe/H] < —2.2 have higher [C/Fe] corr on average than 
stars with [Fe/H] > —2.2. The sample of stars in Fornax 
does not span a large enough range of [Fe/H] to draw 
definitive conclusions on the metallicity dependence of 
[C/Fe] corr . 

We considered the possibility that [C/Fe] corr appears 
to decline with [Fe/H] because our spectra lack the 
S/N or resolution to measure carbon abundances at low 
metallicities. We calculated upper limits for all of the 
stars where we could not measure carbon abundances. 
If the [C/Fe] corr trend appeared because of the loss of 
measurability at low [Fe/H], then the upper limits would 
trace the lower envelope of [C/Fe] corr measurements. In¬ 
stead, the upper limits are distributed fairly evenly in 
[C/Fe] corr . We conclude that the decline is real. 

The decline of [C/Fe] corr with [Fe/H] can be explained 
by the varying contribution of different nucleosynthetic 
sources over time. Figure [TQ] shows theoretical nucle¬ 
osynthetic yields superimposed on the same data as in 
Figure [9] We adopted the AGB yields oflKarakasI (|2010fl . 
the Type II SN yiel ds of INomoto et al.l (120061) . and the 
Type la SN yields of llwamoto et all 1 19991 model W7). 

Type II SNe are the first sources to enrich a galaxy with 
C and Fe. Figure [10] show s the Type II SN yields aver¬ 
aged over a ISalneteiT(119551 1 initial mass function (IMF). 
In pra ctice, th e most massiv e SNe explode first. Accord¬ 
ing to INomoto et al.l (|2006H . a 40 M 0 SN has a C/Fe 
yield 5 times that of a 13 M 0 SN ([C/Fe] = +0.36 ver¬ 
sus —0.36). As a result, the very earliest stars in the 
galaxy may have been able to incorporate the ejecta of 


only the most massive, shortest lived SNe, which have 
the highest [C/Fe] of the Type II SNe. The two stars 
in our sample with [Fe/H] < —3, both in Draco, have 
[C/Fe]corr = —0.39 ± 0.37 and +0.76 ± 0.13. The differ¬ 
ence between these values could indicate separate enrich¬ 
ment from Type II SNe of different masses, or it could 
indicate different levels of AGB pollution. 

The IMF-averaged Type II SN yields are lower than 
the observed [C/Fe] in the majority of dSph stars at 
[Fe/H] < —1.7. These stars require a more carbon-rich 
nucleosynthetic source, such as AGB stars. AGB stars 
were likely polluting the int erste llar medium even at the 
very lowest metallicities. iKarakasI (|2010f ) did not calcu¬ 
late AGB yields at [Fe/H] < —2.3, but AGB stars surely 
existed at lower metallicities. If the high [C/Fe] yields 
persist at [Fe/H] < —2.3, then AGB ejecta could easily 
explain the high [C/Fe] values observed at low [Fe/H] 
in the dSphs. Their ejecta would mix with Type II 
SN ejecta to produce the [C/Fe] ratio observed in low- 
metallicity stars. Because AGB stars have much larger 
[C/Fe] ratios than Type II SNe, it would not take very 
many AGB stars compared to Type II SNe to explain 
the carbon measurements in low-metallicity dSph stars. 

The AGB yields are highly metallicity-dependent. As 
the galaxy’s metallicity increases, the AGB stars eject 
more and more Fe but about the same amount of C. 
From log (Z/Zq) = —2.3 to —0.7, the C output decreases 
by a factor of 2 while the Fe output increases by a fac¬ 
tor of 37. Hence, the metallicity dependence of AGB 
yields could explain some of the decline in [C/Fe] with 
increasing [Fe/H]. 

Type la SNe are the last nucleosynthetic sources to 
turn on in dSphs because they are the end points of some 
long-lived, low-mass stars. The Type la [C/Fe] yield is 
— 1.6, significantly below the yield of Type II SNe and 
AGB stars. As the galaxy evolves, the ratio of Type la 
to Type II SNe increases. Because the abundance ratios 
of the two SN types are drastically different, the trend of 
abundance ratios with [Fe/H] reflects the changing con¬ 
tribution of the two SN types to the g alaxy’s chemical 
evolution (e.g., iGilmore fc Wvsel 11991 1. The elements 
typically used to diagnose this ratio are a elements (O, 
Mg, and Si among others) and Fe. As Figure [12] demon¬ 
strates, C/Fe is also sensitive to the Type II-to-Type la 
SN ratio. Although the heavy production of C in AGB 
stars complicates the interpretation of the C/Fe ratio, 
some of the decline in [C/Fe] with increasing [Fe/H] in 
dSphs must be due to the increasing frequency of Type la 
SNe with time. 

The four dSphs show different [C/Fe] trends. The 
[C/Fe] ratio declines most steeply in Ursa Minor and 
Draco. Under the presumption that an increasing preva¬ 
lence of Type la SNe drives the decline in [C/Fe], 
then these two dSphs would have had a very inefficient 
SFH. The galaxy would have taken about 100 Myr— 
the approximate delay time for a Type la SN—to reach 
[Fe/H] ss —2.5, the metallicity at which [C/Fe] begins to 
decline. Unfortunately, Ursa Minor, Draco, and Sculp¬ 
tor are too old to permit fin e ti me resol u tion in mea¬ 
surements of their SFHs (e.g.. lOrban et al.l[2008 ). so our 
chemical interpretation can not be cross-checked against 
photometrically me asured S FHs. Finally, Fornax has a 
complex SFH (e.g.. Ide Boer et ahll2012H . Its pattern of 
[C/Fe] is unlike the other three dSphs. Some of the most 
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metal-rich stars have high values of [C/Fe] corr . These 
measurements also have some the largest uncertainties, 
largely because the spectrum at 4300 A is blanketed with 
atomic and molecular absorption lines. If these measure¬ 
ments are trustworthy, then we speculate that AGB stars 
are responsible for the high carbon abundances, possibly 
through mass transfer. 

We conclude that the primary causes of the decline of 
[C/Fe] corr with [Fe/H] are (1) the metallicity dependence 
of AGB yields and (2) the increasing ratio of Type la to 
Type II SNe with time. A quantitative interpretation 
of Figures [9] and fTTTl would require a chemical evolution 
model. IKirbv et al.l (1201 ll) developed such a model to 
explain the decline of [ct/Fe] with [Fe/H] in dSphs, in¬ 
cluding the four dSphs discussed here. One prospect for 
future work is to incorporate carbon into the model, pay¬ 
ing close attention to the yields of AGB stars. 

Metallicity does not co mpletel y exp lain the observed 
scatter in [C/Fe]. IGratton et al.l (120001) observed a large 
scatter of [C/Fe] on the upper RGB in low-metallicity 
held stars in the MW halo. However, they did not have 
a large enough sample of unevolved stars at low metal¬ 
licity to determine whether the scatter was imparted to 
the stars at birth or if it was caused by star-to-star vari¬ 
ations in the efficien cy o f car bon astration on the up¬ 
per RGB. ISpite et ahl (j2005T) found an approximately 
constant value of [C/Fe] in unevolved held stars (ignor¬ 
ing C stars), even at low metallicities. This suggests 
that low-metallicity, star-forming gas in the ear l y MW 
halo had a uniform [C/Fe] ratio. However, ISpite et~ahl 
(|2006f) found that evolved stars on the upper RGB did 
not all experience the same level of mixing. In partic¬ 
ular, the [(C+N)/Fe] ratio was larger for evolved stars 
at [Fe/H] < —2.5 than for higher-metallicity stars. In 
CNO-processed material that did not undergo much of 
the high-temperature ON cyc le, t he ] (C+N )/Fe] ratio 
should be constant. Hence. ISpite et all (1200611 uncovered 
a metallicity depe ndence i n the effic iency of mixing. 

In principle, IPlacco et all ’s (120141) metallicity- 
dependent carbon corrections should acco unt f or t he 
variat ion in mixing efficiencies found by ISpite et all 
(|2006l ). However, we still observe a dispersion in 


C/Fejcorr beyond the measurement uncertainty at hxed 
Fe/H . We fit a quadratic relation between [C/Fe] corr 
and [Fe/H] to stars with [C/Fe] corr < +0.7 (to exclude 
stars that may have acquired carbon through binary 
mass transfer). In a manner similar to calculating the 
intrinsic dispersion of [C/Fe] in GCs (Section [4]), we 
calculated the residuals, A, about the quadratic fit 
and divided them by the measurement uncertainties: 
A = ([C/Fe] ?orr - [C/Fe] corr (fit))/<5[C/Fe]. If measure¬ 
ment uncertainties explain the scatter, then the standard 
deviation of A should be 1. Instead, we found standard 
deviations of 1.5 for Sculptor, 1.2 for Fornax, 1.5 for 
Ursa Minor, and 1.6 for Draco. Therefore, all of the 
dSphs show a small residual beyond the measurement 
uncertainties that can not be explained by the decline 
of [C/Fe] with increasing stellar luminosity (Figure [8]) 
or with increasing metallicity (Figure [9]). 

One possible sourc e of the residual scatter is inhomo¬ 
geneous mixing (e.g., iRevaz fe Jablonkall2012l ). If differ¬ 
ent contemporaneous star formation regions formed out 
of gas that was not mixed across regions, then stars at 


the same metallicity and luminosity could have different 
abundance ratios. Also, the calculation of A depends 
on accurate estimates of measurement uncertainty. If we 
have underestimated 6[C/Fe], then we h ave also over ¬ 
estim ated A. Another possibility is that IPlacco et all s 
(120141 ) corrections are imperfect. Errors on the order of 
~ 0.15 dex in the corrections could cause the scatter in 
[C/Fe] corr that we observe. 

5.1. Carbon Stars 

As demonstrated in M15 (Section [4]), stars need not be 
born with their present [C/Fe] ratio. They can acquire 
large amounts of carbon through mass transfer from an 
AGB star. These stars can exhibit [C/Fe] +1. We 
identified such stars in our dSph sample. We found 2 C 
stars in Fornax, 4 in Ursa Minor, 5 in Draco, and none 
in Sculptor. They are represented as red, five-pointed 
symbols in Figure [TJ Most of these stars were discovered 
in previous surveys 0 Table [5] lists all of the C stars 
in our sample along with the the reference(s) and iden¬ 
tification name(s) for their previous discoveries, where 
applicable. They all have extremely strong CN bands, 
as shown in Figure [2] Most of the C stars also have very 
strong C 2 bands visible in the blue spectra. However, 
two of the three Ursa Minor C stars have weak C 2 , and 
C 2 is not visible at all in star Bel80022. It is likely too 
metal-poor and perhaps too warm for C 2 to appear in 
our low-resolution spectrum. 

The resolution of our DEIMOS spectra does not per¬ 
mit a detailed chemical abundance analysis of these car¬ 
bon stars. The CN, C 2 , and CH molecular bands blan¬ 
ket the spectra too heavily. For this reason, their car¬ 
bon abundances are represented as lower limits in Fig¬ 
ure [5] They are not represented at all in Figure [HI be¬ 
cause we could not measure [Fe/H]. Higher resolution 
spectra w ould allo w a detailed chemical analysis (e.g., 
iKarinkuzhi fc Goswamilf2014 ). Almost all of the carbon 
stars we found are in fact bright enough to be observed 
with a high-resolution spectrograph on a large telescope, 
like Keck/HIRES. We reserve this task for future work. 

Figure [3 shows that five of the 11 C stars are redward 
of the RGB. The tendency for carbon stars (or more 
specifically, barium stars, which almost always have ex¬ 
cess carbon) to a ppear red in optical colors is called the 
Bond-Neff effect dBond fc Nefijfl969 ). The source of the 
effect is a depression of blue flux due to molecular blan¬ 
keting of the spectrum. In fact, all five of the C stars 
redder than the RGB show strong C 2 , CH, and CN. 

Our spectroscopic target selection was focused on se¬ 
lecting stars from the RGB. Hence, we may have 
missed many C stars. Although we discuss the C stars 
in our sample, the sample is definitely not complete. 
Other surveys that specific ally searched for C stars (e.g., 
iDemers fc Battinellil 120021) are better able to quantify 
their frequency. Our sample selection for Scul ptor was es¬ 
pecial ly biased against finding C stars. iSkuladottir et alJ 
(12015! ) found the first CEMP star in Sculptor, but it lies 
outside of the area we surveyed. 

The brightest C star in Ursa Minor, Bell0023, is ac¬ 
tually blueward of the RGB. The molecular absorption 


20 IShetrone et alJ (120011 ) identified star 622253 (called 68 in their 
paper) in Draco as a C star. We do not find it to be carbon-rich. 
We measured its carbon enhancement to be [C/Fe] = —0.52±0.11. 
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Table 5 

Carbon Stars 


DSph Name 


Reference 


Reference ID 


Fornax 
Fornax 
Ursa Minor 

Ursa Minor 
Ursa Minor 
Ursa Minor 


Draco 

Draco 

Draco 

Draco 

Draco 


50465 

98788 

Bel60031 

Bel60230 

Bel80022 

Bell0023 


670092 

569394 

607050 

571725 

682522 


this paper 
this paper 

Aaronson et al. (^83) 
Armandroffet^L (1995 
this paper 

^hetroneetah (20(R) 
Canterna &; Schommer 

Zinn (1981) 

Aaronson et al. 

Armandroff et *al. 



Shetroneet al. 

Wi nnick (2003) 
Armandroff et al. (1995) 
Aaronson ct al. (1982) 
Armandro ff et al. (1995) 
Aaronson ct al. (1982) 
Armandroff ct al. (1995) 
Aaronsonet al. (1982) 
Armandroff et al. (1995) 
Azzopardi et al. (1986) 
Armandroff et al. (1995) 


34227 

VA335 

UM 1545 

COS215 

K 

K 

K 

K 

COS215 

461 

3203 

Cl 

J 

C2 

3237 

C3 

578 

C4 


apparently is not strong enough to redden the optical col¬ 
ors. Indeed, the C 2 bands in this star are visible but very 
weak. The star is probably an AGB star, just like 47866 
in M15. It possibly manufactured its carbon by itself 
rather than acquiring it through binary mass transfer, 
though the star is probably too low-mass to dredge up 
carbon. Unfortunately, our spectra do not provide any 
clues as to whether the C stars generated their carbon 
themselves or acquired it externally. 

6. COMPARISON TO THE MILKY WAY STELLAR 

HALO 

The MW stellar halo most likely grows by the minor 
merging of satellite galaxies. As a result, the chemical 
abundance distribution of the halo is the superposition 
of the abundance distributions of all of its constituents. 
The abundances of these individual satellites depend on 
their stellar masses and SFHs. Therefore, the mass dis¬ 
tribution and SFH of the halo’s dissolved satellites can 
be inf erred from the chemical abundances of the halo (see 
iLee et akll2014aD . 

One straightforward way to compare the carbon abun¬ 
dance distributions of the surviving dSphs to the MW 
halo is to measure the frequency of CEMP stars. The fre¬ 
quency depends on metallicity. Sp ecific ally, it decreases 
toward higher metallicity dBeers et all ll99l) . In Sec¬ 
tion 0 we showed that [C/Fe] does indeed decrease with 
increasing [ Fe/H]. After accounting fo r this decrease with 
metallicity, iStarkenburer et ah (l2013i ) presented evidence 
that the CEMP fraction in Sculptor is lower than in the 
MW halo with ~ 90% confidence. This result would sug¬ 
gest that Sculptor had a different chemical evolution or 
perhaps a different frequency of mass-transfer binaries 
than halo progenitors. 

Very few of the stars in our sample would qualify as 
CEMP under the various definitions for “CEMP.” One 
proposed definition is a simple cut in carbon abun dances, 
such as [C/Fe] > +1.0 (|Beers fe Christliebll2005D . How¬ 
ever, this definition disregards the destruction of car¬ 
bon on the upper RGB, which is the location for all of 


the stars in our sample. A more nuan ced definition ac¬ 
counts for the luminosity of the sta r (lAoki et al.ll2007t 
lYong et all 120131 : IPlacco et al.ll2014D . We consider the 
CEMP f ractio n with two different definitions. First, we 
consider lAoki et gill ’s (i2007t ) definition: 

C/Fe] > +0.7 if log(L/L 0 ) < 2.3 

C/Fe] > +3.0 - log(L/L 0 ) if log(L/L 0 ) > 2.3 (3) 

The dotted line in Figure [8] shows the dividing line of 
Equation [3] 

We identified only 6 CEMP stars out of 398 stars 
in our dSph sample using lAoki et all ’s (2007)) defini¬ 
tion of CEMP. Of these, 1 has [Fe/H] < —3.0, 2 have 
[Fe/H] < —2.5, and 5 have [Fe/H] < —1.0. In these same 
metallicity ranges, the CEMP fractions are (50 ± 50)%, 
(7.1 ± 5.1)%, and (1.3 ± 0.6)%. However, our sample is 
heavily biased against extremely carbon-rich stars. First, 
we selected against stars redward of the RGB, where 
C stars are generally found. Second, our technique of 
measuring [Fe/H] in the range 6300-9100 A fails for ex¬ 
tremely carbon-rich stars, in which CN blankets much 
of that spectral region. We found 11 C stars, but we 
did not measure their metallicities. As a result, they are 
not included in the above estimate of CEMP frequency. 
Therefore, the frequencies of carbon stars we found in 
dSphs should be regarded strictly as lower limits. 

We also counted CEMP stars using [C/ Fe]„„^, w hich 
is cor rected for luminosity according to IPlacco et all 
(|2014T> . One of IPlacco et all ’s proposed definitions for 
CEMP is [C/Fe] corr > +0.7 (dotted line in Figure [p]). 
The cut does not need any dependence on luminosity 
because [C/Fe] corr already depends on luminosity. Un¬ 
der this definition, we found only 3 out of 398 stars. 
They have [Fe/H] = -3.22 ± 0.20, -2.51 ± 0.26, and 
-2.07 + 0.11. 

The CEMP fraction of the halo is significantly higher 
than the lower limits we determined for dSphs ex¬ 
cept in the [Fe/H] < —3.0 range, where we ob- 
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Figure 11. The trend of [C/Fe] corr versus [Fe/H] in dSphs com¬ 
pared to stars in the MW (top panel: Gratton et al. 20(50 and 
I Spite et al.l 120051 : bottom panel: ILee et aD I2013f ) . In the bot¬ 
tom panel, the grayscale for MW stars is logarithmic in number 
density. The bottom panel shows only main sequence turn-off 
stars, which were selected in the same manner as in ILee et aLl : 
5600 K < T e ff < 6700 K. These stars have not altered their sur¬ 
face compositions, so their [C/Fe] measurements do not need to be 
corrected. 


served only 2 stars. In the halo, the CEMP fraction 
is approximat ely 13% at [Fe/H] < —2 and 25% at 
[Fe/H] < —3 dCo hen et all 120051 i Marsteller et all 120051 


Lucatello et alT 20061 iFrebel et all 120061 iCarollo et al 

20121 lYong et al.l 120131: ILee et all 120131 iPlacco et al 


2014f) . Many of the halo’s CEMP stars have [C/Fe] > +2. 
These stars would have been excluded from our sample. 

Nonetheless, we can still compare the general trend 
of [C/Fe] versus [Fe/H] for stars with [C/Fe] < +1 be¬ 
tween dSphs and the halo. Although our sample misses 
many of the stars with [C/Fe] > +1, the remaining stars 
show a trend different from the MW halo, as shown in 
Figure fill The top panel compares our measurements of 
[C/Fe] corr in dSphs com pared to high-resolution measure- 


ments in the MW halo dGratton et al. 


2005). The sample of lOratton et al 


2000t iSoite et all 


includes dwarfs 


and giants, and t he sample of ISoite et al.l is only gi¬ 
ants. We applied IPlacco et ahl ’s (120141) corrections to 
the [C/Fe] measurements in both samples. The bottom 
panel compares our measurements to SEGUE main se¬ 


quence turn-off stars (|Lee et alJl2013l l. measured by the 
Sloan Extension f or G a lactic Unde rstanding and Explo¬ 
ration (SEGUE. lYannv et al.ll2009li . We plot the turn-off 
stars because their carbon abundances have not yet been 
affected by stellar astration. We adjusted all of the MW 
[C/Fe] measurements to match the solar abundance ra¬ 
tios we adopted (see Section l5.4.3l) . Ignoring the CEMP 
stars, the [C/Fe] pattern of the general halo population 
is flat at [C/Fe] ~ +0.1 until [Fe/H] = —0.7, whereupon 
it drops to [C/Fe] ~ —0.2 at [Fe/H] = 0. In contrast, 
the dSph stars have already dropped to [C/Fe] ~ —0.4 
by a metallicity of [Fe/H] = —2. Note that the floor 
in the SEGUE data at [C/F e] ~ — 0.3 might be artificial. 
The lowest value of [C/Fe] in lLee etldl ’s grid of synthetic 
spectra is —0.5, and the minimum upper limit values (not 
plotted in Figure fTlTi are mostly about 0.2 dex away from 
the grid limit. Regardless, the SEGUE s amp l e is q ual¬ 
itative ly c onsistent wi t h the samples of iGratton et al.1 
(1200011 and [Spite et al.1 (l2005lh 
ILee et ah ( 2014bfl interpreted the halo’s CEMP frac¬ 
tion in the context of binary mass transfer. They found 
that carbon donation from AGB stars could explain the 
present [C/Fe] distribution for giant stars. However, 
that mechanism is primarily responsible for large car¬ 
bon enhancements ([C/Fe] > +1). In Section [5] and 
Figure [TU] we interpreted the trend among non-CEMP 
stars in the context of galactic chemical evolution. The 
different [C/Fe] slopes between the dSph and halo sam¬ 
ples is reminiscent of the different slop es in [a/Fe] (e.g., 
iShet.rone et, al.lll998al: I Venn et, al.ll2004H . The [a/Fe] dis¬ 
tribution of the halo shows a downturn (a “knee”) at 
[Fe/H] = —0.7, a feature that is also visible in [C/Fe] 
(Figure |TT|). The knee corresponds to the metallicity 
at which Type la SNe began to pollute the MW. Be¬ 
cause Type la SNe have low yields of both [a/Fe] and 
[C/Fe], their prevalence beginning at [Fe/H] = —0.7 can 
explain the abundance pattern for both abundance ra¬ 
tios in the halo. However, the knee happens at lower 
metallicities for dSphs. The decline of [C/Fe] with in¬ 
creasing [Fe/H] begins at [Fe/H] < —2.5 for Ursa Mi¬ 
nor and Draco. The lower metallicity of their [C/Fe] 
knees indicates that Type la SNe were active at much 
lower metallicities in dSphs than in the halo. The same 
trend is seen in dSphs ’ [a/Fe] ratios (e.g., lYenn et all 
12004 iKirbv et al.ll201 lT) . The halo reached a much higher 
metallicity than the dSphs before Type la SNe turned on. 
Furthermore, both [C/Fe] and [a/Fe] in dSphs decline to 
values well below the solar ratio, indicating a heavy dose 
of Type la SNe relative to the halo, which has larger 
values of [C/Fe] and [a/Fe]. Our [C/Fe] measurements 
are further proof that the chemical evolution of surviving 
dSphs was much slower than the progenitors of the MW 
halo. 


7. SUMMARY 

Carbon is a complex element. It has diverse nucle- 
osynthetic sources. AGB stars produce copious amounts 
of carbon, but their C/Fe ratios depend sensitively on 
metallicity and stellar mass. Type II SNe produce some 
carbon, whereas Type la SNe produce very little car¬ 
bon but a great deal of iron. The timescale for all three 
sources is different, which makes the C/Fe ratio change 
as the galaxy evolves in a way that depends on the details 
of the SFH. Mass transfer from past AGB stars further 
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complicates the interpretation of the carbon abundances 
of present-day stars. In particular, some stars in mass- 
transfer binaries can have orders of magnitude more car¬ 
bon than single stars of the same metallicity. In addition 
to multiple sources, carbon also has a sink. Red giants 
astrate carbon on the upper RGB. The rate of astration 
depends on the metallicity and luminosity of the star. 

In order to explore the evolution of carbon, we ob¬ 
served 745 red giants with the Keck/DEIMOS spectro¬ 
graph. We generated a grid of model spectra based on a 
line list that we checked against spectra of the Sun and 
Arcturus. We measured carbon abundances by fitting 
the model spectra to the G band at 4300 A. We were 
able to measure [C/Fe] for 552 stars. We validated our 
measurements by comparing a subset of our observations 
to published measurements of carbon abundances for the 
same stars. We found offsets between our measurements 
of [C/Fe] and the literature measurements. However, the 
offset depended on the source of the measurements. Our 
measurements agreed most closely with those based on 
atomic and molecular line lists most similar to our own. 

Globular clusters provide controlled environments to 
study the evolution of carbon because they generally 
have single ages and metallicities. We measured [C/Fe] 
in 154 red giants in the GCs NGC 2419, M68, and M15. 
We observed the well-known decline of [C/Fe] on the up¬ 
per RGB starting at the luminosity function bump. We 
found that carbon does not partici pate in the u nusual 
abundance patter ns in NGC 2419 (|Cohen fe Kirbvll2012t 
iMucciarelli et al .11 20 121) . We also reported the discovery 
of a carbon-enhanced AGB star (47866) in M15. 

We also explored the evolution of carbon in MW dSphs 
by examining 398 carbon abundances and 182 upper lim¬ 
its for stars confirmed to be members of Sculptor, Fornax, 
Ursa Minor, and Draco. As expected, the [C/Fe] distri¬ 
bution is more complex in dSphs than in GCs. Unlike 
GCs, dSphs show a dispersion in [C/Fe] at fixed stellar 
luminosity. The dispersion depends on [Fe/H]. Specif¬ 
ically, stars with lower [Fe/H] have higher [C/Fe] than 
stars at the same luminosity with higher [Fe/H], We in¬ 
terpreted the trend with metallicity as an evolution in the 
influence of the three sources of carbon. First, Type II 
SNe seeded the dSph with the first C and Fe. Very soon 
after, the first AGB stars produced a great deal of C rel¬ 
ative to Fe. Finally, Type la SNe produced very little C 
but a large amount of Fe after ~ 100 Myr. The result is 
a decline in [C/Fe] with increasing [Fe/H] that roughly 
tracks t he decl ine in [a/Fe], also due to Type la SNe 
(e.g., iKirbv et al.ll2011[) . 

Our color selection biased the sample against very 
carbon-rich stars. Nonetheless, we found 11 C stars (8 
previously known) in three dSphs. These stars are either 
AGB stars that are producing carbon now or the results 
of binary mass transfer from past AGB stars. Because of 
our selection bias, we calculated only lower limits to the 
CEMP fraction in dSphs. These lower limits are mostly 
below the CEMP fraction for the MW halo. Therefore, 
we have found no tension for the theory of hierarchical 
assembly, wherein the halo builds mass t hrough minor 
merg ing of satellites like dSphs (e.g.- lBullock fc Johnstonl 
120051) . However, this question would be addressed better 
by a spectroscopic survey of dSphs with a color cut more 
inclusive of stars redder than the RGB. 


Our sample is still useful for comparing the trend of 
[C/Fe] versus [Fe/H] for carbon-normal stars. We found 
that the knee in the [C/Fe] distribution occurs at a l ower 
metallicity in dSphs than in the MW halo dLee et all 
120131 1. This difference in metallicity is reminiscent of the 
different [a/Fel knees between dSphs and the MW halo 
(iVenn et all 120041) . The similarity in behavior between 
[C/Fe] and [a/Fe] supports our hypothesis that Type la 
SNe are responsible for the decline in [C/Fe]. The [a/Fe] 
and [C/Fe] distributions together indicate that the pro¬ 
genitors of the MW halo had a more efficient SFH, which 
reached a higher metallicity before Type la SNe began 
to produce Fe in earnest. 

Our sample of 398 measurements doubles the num¬ 
ber of carbon abundances known in dSphs. Each dSph 
has enough carbon measurements (181 in Sculptor, 45 
in Fornax, 41 in Ursa Minor, and 131 in Draco) for 
analysis with a numerical chemical e volution model (e.g., 
iRevaz et al.ll2009t IKirbv et al.ll2011lf . Although we have 
presented a qualitative interpretation for the observed 
carbon abundance trends, a numerical model will give a 
quantitative timescale for star formation. This chemi¬ 
cal timescale can be compared to photometric measure¬ 
ments of SFH. Further more , nume rical models can possi¬ 
bly c onstrain the IMF dLee et al.ll2014bl) or binary frac¬ 
tion dLucatello et all I2005D in dSphs. We invite inter¬ 
ested modelers to make quantitative models of the car¬ 
bon abundances in Table []] to reveal the history of star 
formation in these four dSphs. 
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Table 4 

Carbon Abundances 


System 

Name 

RA 

(J2000) 

Dec 

(J2000) 

log L 
(£©) 

S/N 

(A- 1 ) 

TeS 

(K) 

log 9 

(cm s -2 ) 

[Fe/H] 

(dex) 

[a/Fe] 

(dex) 

[C/Fe] 

(dex) 

[C/Fe] corr - 
(dex) 

S 2 (CH) 

(mag) 

C star? 

Ursa Minor 
Ursa Minor 

Bel60031 

Bel60021 

15 07 38.73 
15 07 50.88 

+67 13 56.4 
+67 11 12.9 

2.94 

2.63 

99 

4495 

1.30 

-2.43 + 0.11 

+0.08 + 0.15 

-0.65 + 0.16 

-0.01 


c 

Ursa Minor 

Bel60469 

15 07 55.52 

+67 11 30.2 

1.74 

32 

4772 

2.33 

-2.78 + 0.25 


< +0.44 

< +0.45 

1.75 + 0.33 


Ursa Minor 
Ursa Minor 

Bel60230 

Bel60123 

15 07 59.77 
15 07 59.77 

+67 15 28.8 
+67 11 41.4 

2.28 

2.10 

56 

4826 

1.96 

-2.36 + 0.13 

+0.35 + 0.19 

< +0.06 

< +0.13 

1.60 + 0.11 

c 

Ursa Minor 

Bel70584 

15 08 08.20 

+67 06 43.8 

1.62 

24 

4899 

2.49 

-2.30 + 0.21 

-0.19 + 0.68 

< +0.55 

< +0.56 



Ursa Minor 

Bel70182 

15 08 08.29 

+67 04 41.0 

2.00 

25 

4803 

2.06 

-2.02 + 0.16 

-0.17 + 0.46 

< +0.90 

< +0.94 



Ursa Minor 

Bel60017 

15 08 09.06 

+67 09 21.6 

2.63 

70 

4526 

1.33 

-2.07 + 0.11 

+0.17 + 0.13 

+0.31 + 0.27 

+0.73 

2.13 + 0.05 


Ursa Minor 

Bel50017 

15 08 10.51 

+67 17 07.1 

2.60 

111 

4529 

1.35 

-2.18 + 0.11 

+0.04 ± 0.08 

-0.70 + 0.21 

-0.10 

1.59 + 0.05 


Ursa Minor 

Bel60819 

15 08 12.07 

+67 10 15.5 

1.27 

13 

5198 

2.93 

-1.85 + 0.33 


< -0.26 

< -0.25 

1.54 + 0.17 



References. — Photomet ry comes from the following sour ces: NGC 2419 and M15 H Stetson! 20( Kj) , M68 jXlcaino et alJ 119901 : IWalkeH 11994) : IStets on 2000), Sculptor 
d Westfall et al. I [20061) . Fornax (IStetson et al.lfl998i ). Ursa Minor dBellazzini et all 1 20021b Draco fSega 1 et alJI2007T) . 

Note. — (This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. The 
online version of the table has additional columns for photometric magnitudes.) 
a [C/Fe] corrected for stellar evolutionary state according to lPlacco et alJ (|2014l) . 
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